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Recently,  the  Army  has  shewn  interest  in  minimizing  "collateral 
damage*"  to  civilian  areas  and  to  friendly  troops.  This  interest  has 
created  a need  for  facilities  to  permit  "full  scale"  testing  under 
simulated  low-level  blast  loadings  caused  by  explosive  charges.  ("lull 
scale"  walls  would  he  of  realistic  sizes  wherein  regular  construction 
materials  and  methods  are  useable.)  The  need  for  such  facilities 
apparently  has  not  been  met  satisfactorily.  facilities  exist,  e.g. 
the  DASACON  conical  tube  at  Dahlgrcn,  Virginia,  permitting  full  scale 
testing;  yet  operationally  they  may  have  shortcomings:  expense  of 
operation  in  terms  of  manpower  and/or  energy  requirements  or,  in- 
adequate simulation  of  the  waveform  characteristics  of  conventional 
explosive  blast  waves.  In  response  to  an  Army  request  involving 
such  full  scale  testing,  we  have  studied  a possible  modification  of 
the  2.4m  (81)  tube  portion  of  the  existing  BRL  Dual  Shock  Tube 
Facility1’2,  Figure  1 from  Reference  1,  to  give  the  desired  testing 
parameters  - acceptance  of  walls  to  3.7m  (12')  height,  subjected  to 
peaked  pressure  pulses  up  to  70  kPa  (10  psi)  of  'v  50ms  positive 
duration. 

Testing  beyond  the  open  end  was  not  considered  because  of  the 
jet-like  characteristics  of  the  flow  emitted.  Close  in,  pressures 
are  uneven  across  a planar  surface  immersed  normal  to  the  flow,  a 
result  predicted  from  a simple  model,  and  found  in  early  experiments 
by  Bertrand2.  Farther  away,  when  the  shock  is  more  planar,  the  desired 
pressure  levels  are  not  assured  under  the  tube's  operating  limits. 

The  modification  considered  is  basically  an  adjustment  on  the 
driver/driven  tube  lengths  and  test  station  location  to  give  the 
desired  pressure  wave  form,  and  is  a somewhat  natural  extension  of 
the  classical  shock  tube  operation  cycle,  but  with  an  area  change. 

The  peaked  pressure  waveform  is  obtained  by  testing  at  stations  within 
a distance  interval  defined  by  the  intersection  of  the  leading  and 
trailing  waves  of  the  breech-reflected  expansion  fan  with  the  primary 
shock.  The  interaction  of  the  expansion  waves  with  the  shock  weakens 
it,  thereby  giving  the  desired  decaying  pressure  pulse.  The  flow  in 
an  equivalent  straight  shock  tube  is  used  to  follow  the  wave  pro- 
cesses. Specifying  the  desired  positive  duration  of  the  peaked  wave 
then  determines  the  appropriate  tube  lengths,  for  the  straight  tube, 
then,  by  imposing  a simple  physical  requirement,  these  tube  lengths 
can  be  carried  over  to  determine  the  desired  area-changed  tube  lengths. 


1.  B.  P.  Bertrand,  "BRL  Dual  Shook  Tube  Facility, " BRL  MR  2001, 
U.  S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving 
Ground,  MD,  August  1969.  AD  693264. 
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The  driving  pressures  for  the  area-changed  shock  tube  are  computed 
following  the  methods  of  Alpher  and  White,3  who  also  treated  area- 
changed  shock  tubes  but  whose  interests  were  in  stronger  shocks,  with 
the  driver  tube  looking  into  a smaller  area  driven  tube.  The  principal 
features  of  our  calculations  apparently  have  not  previously  been  done 
explicitly  for  the  subsonic  flows  occurring  here  and  up  to  the  wave 
interaction  regions  involved.  Results  are  presented  of  calculations 
for  the  driver  pressures  for  the  2.4m  tube  for  several  test  condi- 
tions, involving  two  different  (larger)  areas  of  driven  tube. 


II.  DETERMINATION  OF  DRIVER  PRESSURE 

The  flow  situation  for  the  area-changed  tube  is  depicted  in 
Figure  2,  along  with  its  wave  diagram  or  x-t  plot.  Following  Alpher 
and  White,  we  may  write  the  expression  for  the  diaphragm  pressure 
ratio  after  the  diaphragm  is  broken  and  flow  processes  are 

established,  as  (ideal,  one-dimensional  flow  throughout) 

P4c  _ P4c  p3a  £3  £2^ 

P1  " P3a  P3  P2  P1 

The  pressure  ratios  in  the  expression  are  interpreted  to  be 

P4c^P3a  " t^ie  Pressure  (ratio)  required  to  expand  unsteadily 
the  gas  in  region  4 from  rest  to  a Mach  number 


P3a/p3 


P3/p2 


- the  pressure  required  to  expand  steadily  the  gas 

from  M,  to  M,*: 

3a  3 

- the  pressure  across  the  contact  surface  separating 
driver-  from  driven-gas,  with  p,  = p_; 


P2/Pl 


- the  pressure  required  to  compress  the  test  gas  in 
region  1 through  the  shock  of  strength  M^  or  Ms> 


Assuming  the  processes  are  isentropic  and  accounting  for  the  steady  or 
unsteady  nature,  we  may  write  Equation  (1)  as 


3.  R.  A.  Alpher  and  D.  R.  White , "Flow  in  Shook  Tubes  with  Area 
Change  at  the  Diaphragm  Section, " J.  Fluid  Meah  3,  457-70 
(1958). 


*For  subsonic  flow  M < 1 , a steady  expansion  is  more  efficient  in  ^ 

the  conversion  of  thermal  to  kinetic  energy  than  an  unsteady  expansion. 

4.  E.  L.  Resler,  Shao-Chi  Lin , and  Arthur  Kantrowitz,  "The  Production 
of  High  Temperature  Gases  in  Shock  Tubes,"  J.  Appl.  Phys.  23^  1390-99 
(1952).  10 
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i,  Y,-l  2+  (Y.-l)  M " , P2 

iii  = (I*  -4-  »3a>  l 5 V‘  f-  l2) 

P1  ' 2*  (V‘)  H3»  * 

In  regular  shock  tube  testing,  the  gases  in  regions  1 and  4 are 
usually  specified,  as  is  the  desired  test  overpressure;  thus  a driver 
pressure  can  be  determined.  In  the  area-changed  tube  the  unknowns 

are  p as  well  as  M_  and  M,  . Equations  to  connect  these  are  the 
1 4c  3 3a 

expression  for  isentropic  nozzle  flow  in  terms  of  the  nozzle  areas. 

1 V1 

2 2 Y4-:1 

A,  _ Mj  2»Cy4-D  M3a 

A1  M3a  2*(y4-1)  Mj2 

and  an  expression  connecting  with  the  driven  gas: 

= U3  _ U3/al  _ U2/al  f4) 

M3  a3  a^/a^  (a^/ a?a)  .(a^/ a^  (a^/'a^J 

where  a similar  decomposition  into  steady-unsteady  processes  as  in 
Equation  (1)  is  made,  and  results  in: 

(u  /a  )(a  /a  ) 

M = — - — — — 

2+(Y4-l)  M3a2  J 2 

[ 1 2±_]  2 = (5) 

2+(y4-1)  M3"  2+ C Y 4 — 1 ) M3a 

We  now  combine  Equation  (5)  with  Equation  (2)  and  rewrite  (2)  as 


— = t 

n L 


(u2/a1)  (a1/a4) 


P2/pl 


Equations  (5),  (6)  and  (3)  may  then  be  used  to  determine  the  required 
driver  pressure  P4c>  knowing  the  desired  testing  pressure  P2/P^ 

the  nozzle  area  ratio  AVAj.  Solution  techniques  - e.g.  the  Newton- 

Raphson  method  or  iteration  - may  be  used  on  this  system  of  equations. 
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We  have  used  a graphical  procedure  offering  some  simplicity  and 
utility  to  the  worker  at  the  site.  This  is  now  illustrated.  From 
liquation  (3)  we  have 


(A.± 


A.  M, 

4 aa 


Vi 

V1 


A1  M3 


t2*  (y1*  1/2 

2+  (Y4-l)  M2 


which,  on  substituting  into  Equation  (5),  gives 

(u2/a1)  (a2/a4) 


(\  M 

\A1  M3  ) 


V1 

Y4+1 


2+CY4"1)  M3a 


(7) 


Or  rearranging 


V1 

V1 


2 cyv 


Y4  + 1 


(u2/a1)  (a1/a4)  3 


2+  ‘V1)  M3a 
Y4-1 


V1 


3 a 


(8) 


The  right  hand  side  of  Equation  (8)  determines  a curve  Z plotted 

J a 

vs  Mach  number;  similarly,  the  left  hand  side  determines  a curve 
(a  straight  line  on  a log-log  plot).  Thus  for  given  p^p^  ~ 
or  (a1/a4)  - and  A4/A^,  the  points  Z^  = Z^&  determine  the 

Mach  numbers  M and  M^  upstream  and  downstream,  respectively,  of 
the  nozzle,  and  in  particular  M^  of  Equation  (6)  and,  thereby,  the 
required  driver  pressure  P4c* 


Some  sample  values  have  been  calculated  for  their  specific 
interest  and  to  illustrate  the  procedure.  An  attractive  feature  is 
that  different  area  ratios  and  test  pressures  are  reflected  only  in 
a change  in  (straight  line)  intercept  (for  the  same  gases)  of  the  Z_ 

curve.  The  area  ratios  considered  were  for  2.4m  (8')  diameter 
driver  tube  and  Case  (a)-3.7m  (12')  and  Case  (b)  - 4.6m  (15') 
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= 73  kl’a  (10.6  psi  ) 


diameter  driven  tubes,  with  i test  overpressure  p } 


liases  had 


Y4  = Y1  = 7/5;  a4  = :‘l 


Case  (c)  is  similar  to  (a)  but  p,  = 34.5 

for  the  various  quantities  entering  into 
are  tabulated  in  Table  IA. 

kPa  (5.0  psi).  The 

the  construction  of 

values 

the  graph 

table: 

IA  - PARAMETERS  FOR 

TEST  FLOWS 

Case 

(A4/Ai) 

U2/al  M1 

2(A4/A1) 

1/^(>(u2/a1) 

P2/P] 

a 

0.444 

0.4022  1.27 

4 

. 346 

1.72 

b 

0.284 

0.4022  1.27 

4 

.033 

1.72 

c 

0.444 

0.219  1.14 

7 

.982 

1.34 

The  graphs  for  these 
Mach  number. 

cases  are  plotted  in 

Figure  3 

as 

functions  of 

For  the  three  cases  the  Mach  numbers 
Table  IB  as  well  as  the  driver  pressures 
from  Equation  (6). 

detcrmined,ajre  tabulated,  in 
(absolute)  as  determined 

TABLE 

IB  - CALCULATED  NOZZLE  MACH  NUMBERS  ANL)  DRIVER 

PRESSURES 

Case 

(kPa) 

P4c 

(psi)  'v 

M_ 

.*>a 

a 

800 

116 

0.6 

0.50 

b 

1564 

227 

0.68 

0.55 

c 

252 

37 

0.40 

0.25 

From  these  sample  cases,  it  is  seen  that  a 2.4m  diameter  driver 
section  coupled  to  a 3.7m  diameter  driven  section  to  give  test  over- 
pressures up  to  70  kPa  (10  psi)  is  well  within  the  1100  kPa  operating 
limits  of  the  present  TBD  2.4m  shock  tube.  A 4.6m  driven  section 
at  this  overpressure  could  not  be  accommodated,  however.  For 
reference  we  have  considered  a range  of  test  over-pressures  p_,^  and 

their  required  driver  piessures  p for  the  2.4m:  3.7m  (8':12') 
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area-changed  tube.  Results  are  plotted  as  Figure  4.*  Maximum  p is 

seen  to  be  'v  93  kPa  (13.5  psi)  for  the  1100  kPa  driver  operating 
1 imit . 

In  addition,  some  cases  have  also  been  treated  involving 
the  1.7m  (5  1/2')  diameter  tube  of  the  BRL  Dual  Shock  Tube  Facility, 
to  illustrate  the  obtainable  test  pressures,  along  with  associated 
driver  pressures.  Wave  durations  were  not  worked  out.  The  1.7m  tube 
lias  provision  for  heating  of  the  driver  gas  to  eliminate  the  contact 
surface  discontinuity,  i.e.  for  tailoring  of  the  interface. 

For  small  disturbances,  identity  of  the  acoustic  impedance  pa 
across  the  interface  assures  that  waves  are  transmitted  without  change. 
For  the  same  gas  on  either  side  of  the  interface,  the  requirement 
reduces  to  identity  of  the  sound  speeds.  Hence,  also,  in  the  tailored- 
interface  mode 

u,/a_,  = u_/a^  = M^. 

With  the  desired  test  condition  (hence  M^)  known,  one  may  enter 

the  plot  of  Figure  3 to  select  a particular  pair  (there  is  a range) 

M„  > M„  for  which  Z = Z , as  required  by  Equation  (8)  in  the  solution 

for  p . Through  this  point  (M  ,Z^)  is  passed  a straight  line 

paralleling  previously  computed  straight-line  Z^  curves  for  other 

conditions;  and  this  gives  the  intercept  value,  from  which  the  driver 
sound  speed  ratio  a. /a,  is  determinable.  Consequently,  the  driver 

gas  temperature  or  the  heating  may  be  determined.  Then  the  driver 
pressure  including  the  heating  is  again  determined  from  Equation  (6) . 
Three  sample  cases  are  given  in  the  following  Table  II,  for  the  1.7m 
tube  coupled  to  a 3.7m  diameter  driven  tube,  with  T = 300  K. 


* Because  of  the  relatively  Iccrge  exponent  occurring  in  Equation  (6), 
small  rounding  errors  in  the  early  stages  of  calculations  can  yield 
variations  from  our  final  results.  The  accuracy  conforms  to  the 
approximate  nature  of  the  inquiry,  and  does  not  appear  to  be  the 
limiting  feature  in  the  tube  modification. 
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A.  Wave  Advance  in  the  Straight  Shock  lube 

We  turn  now  to  questions  on  the  positive  duration  of  the  desired 
peaked  wave.  As  mentioned  in  the  Introduction,  we  rely  on  the  inter- 
action between  the  incident  shock  S and  the  reflected  rarefaction  wave 
R^  from  the  breech  to  weaken  the  shock  pressure,  thus  giving  the  peaked 
pressure  wave  form.  The  wave  diagram  of  Figure  5,  for  a straight  shock 
tube,  illustrates  a possible  set  of  events  and  the  pressures  at  various 
stations  x.  For  sufficiently  strong  rarefaction  waves  the  shock  reduces 
to  one  of  the  characteristics  of  the  rarefaction  wave  exiting  from  the 
interaction  region;  for  weak  rarefaction  waves  overtaking  the  shock,  a 
weakened  shock  with  lowered  shock  pressures  emerges-1.  In  both  instances 
we  may  take  the  peaked  wave's  positive  duration  as  the  interaction  time 
of  the  shock  with  the  reflected  expansion  fan  from  the  driver  breech. 

The  problem,  therefore,  is  to  determine  this  time  from  the  operation 
cycle  of  our  area-changed  shock  tube. 

As  a point  of  departure,  we  make  use  of  calculations  involving  a 
straight  shock  tube  chosen  to  give  the  same  downstream  conditions  as 
the  area-changed  tube.  The  interaction  problem  and  the  relevant  times 
will  be  solved  for  the  straight  tube.  Then  asking  for  a simple  physical 
requirement  enables  a transforming  of  the  computed  numbers  to  the 
actual  tube. 

Most  of  the  features  of  the  wave  interaction  problem  in  the 
straight  shock  tube  are  well  known.  The  interactions  delineate  in  the 
(x,t)  - plane  a number  of  regions  of  shock  tube  flow.  These  have  been 
labeled  in  the  Figure  5 according  to  standard  practice.  We  focus  on 
the  interactions  involving  the  reflected  leading  and  trailing  waves 
of  the  expansion  fan  connecting  straight  shock  tube  regions  4 and  3. 
Where  and  when  these  overtake  the  incident  shock  wave  tell  a)  where  a 
test  wall  should  be  placed  to  experience  a peaked  pressure  pulse  and 
b)  what  positive  duration  is  to  be  expected.  Conversely,  if  the 
positive  duration  is  specified,  the  problem  formulation  allows  driver 
and  driven  lengths  to  be  determined. 

Thus,  considering  the  propagation  of  the  forward-facing,  reflected 
rarefaction  waves  R,  we  note  that  they  will  first  interact  at  the 
contact  surface  K.  This  interaction  may  lead  to  weak  reflected 
rarefaction  or  compression  waves,  but  always  to  a transmitted  rare- 
faction wave.  The  transmitted  rarefaction  wave  thus  continues  onward, 
overtakes  and  interacts  with  the  primary  shock.  This  interaction  gives 
rise  to  four  possibilities  depending  on  relative  strengths  of  the 


5.  I.  I.  Glass  and  J.  G.  Hall,  Handbook  of  Supersonic  Aerodynamics, 
Section  18,  Shock  Tubes , NAVOHD  Report  1488  (Vol.  6)  1959. 
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interacting  waves:  strong  rarefaction  wave  relative  to  shock  gives  Ij 
reflected  rarefaction  or  J)  compression  waves  and  a transmitted 
rarefaction  wave;  or  weak  rarefaction  wave  relative  to  shock  gives  5}  re- 
flected rarefaction  or  4)  compression  waves  and  a transmitted 
(weakened)  shock.  The  different  initial  conditions  required  to 
achieve  these  cases  could  he  a means  of  achieving  different  pressure 
pulses  with  changed  wave  slopes,  if  desirable,  for  this  study  it 
turns  out  that  the  last-named  interactions  did  not  have  to  be  considered 
in  such  detail. 

15.  Determination  of  Wave  Speeds 

The  progress  of  the  R v wave*  from  the  breech  after  reflection 
requires  knowledge  of  its  wave  speed  C through  the  various  uniform 
regions.  The  wave  speed  is  obtained  from  the  Riemann  invariants 
governing  the  simple  flow.  Across  the  backward  facing  R ^ wave 
connecting  regions  3 and  4: 


:4-l  a4  = U3  + Y 4 — 1 “3  " “3  u>,3  T y^l' 


a = a_  (M 


u = 0 
4 


V1 


14,  + 1 


u,  = r 

3 Y„-l 


(a,  - a„) 

v 4 y 


-y , 


V1 


— = (— ) 
P,  a, 


Have  Speeds: 


C 3 = u,-a3>  C 4 = -a^ 


Across  the  forward-facing  R wave  connecting  3 with  5: 


*The  notation  is  R = rarefaction,  S = shook,  and  arrows  indicate 
direction;  and  C 4 = wave  speed  of  backward- facing  wave  in  region  4. 


Wave  speeds:  C+  3 = + a^;  C+5  = + a,. 

{weak  KR_^ 
or 

weak 

compress,  wave,  S^KR^ 

If  we  assume  the  result  R^  KR^,  our  calculations  lead  to  a 

contradiction  in  the  pressures  across  the  wave  as  well  as  in  the  slopes 
of  leading  and  trailing  wave.  In  this  regard,  the  criteria  and  out- 
come for  this  interaction,  as  stated  by  Courant  and  Friedricks,6 
page  180,  give  an  incorrect  result  here,  whereas  those  by  Glass  and  Hall,' 
page  95-6,  and  Landau  and  Lifshitz7  (for  acoustic  waves)  page  255-6, 
predict  the  results  we  indeed  find. 

The  interaction  results,  then,  in  reflected  compression  waves 
coalescing  to  a shock  and  a transmitted  rarefaction  wave.  For  the 
weak  waves  expected  here,  the  compression  waves  have  the  same  properties 
as  weak  shocks.5  Hence  we  assume  the  shock  wave  equations  are 
applicable  and  also  that  entropy  changes  across  the  "forming  shock" 
may  be  neglected,  i.e.,  the  process  is  also  assumed  to  be  isentropic. 

With  these  assumptions,  regions  5 and  6 are  connected  by 


6.  R.  Courant  and  K.  D.  Friedrichs , Supersonic  Flow  and  Shock  Waves. 
Interscience  Publishers , Inc.,  New  York  (1948). 

7.  L.  D.  Landau  and  E.  M.  Lifshitz , Fluid  Mechanics , Pergamon  Press 
Ltd.,  London,  1959. 
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A nJ~ 

V“s  - as  t - 


2/y, 


5 »-4  + 1)  W + fVU 

2V, 


and 


= A)  V1 

P5  a5 


using  the  isentropie  assumption. 


On  the  other  side  of  the  contact  surface,  across  the  transmitted 
forward- facing  R connecting  regions  2 and  7,  one  has 


2 2 
u„  - a„  = u„  - — — r a 


2 Yj-1  2 7 y1~1  7 


and 


2Y, 


V1 


7 “7 

= (r1) 


where  the  known  conditions  for  region  2 are  obtained  from  the  shock 
strength  p^/Pj  or  Ms>  i.e.,  from  tables  of  straight  shock  tube 

properties. 


The  first  equation  is  recast  into  more  useful  form: 


V1 


2Y, 


u_  = u.  + 


ao  [(~) 


7 “2  Yj-1  2 lvp2 


-1]  , 


and  making  use  of  equalities  across  the  contact  surface 


U6  = u7*  P6  = P7:  P3  = P2 
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results  in 


yielding  a single  equation  for  P^/P5  in  terms  of  previously  determined 

quantities.  Solution  then  permits  the  calculation  of  particle  and 
sound  speeds  u^,  a&  and  u^,  a^  applicable  for  these  regions. 

The  transmitted  leading  and  trailing  waves  of  the  expansion  then 
propagate  onward  with  velocities  (u2+a2)  and  (u^a^),  respectively, 

to  interact  with  the  shock  giving  the  locations  within  which  one  may 
expect  a peaked  wave,  and  also  determining  its  duration.  Specific 
cases  of  test  conditions  must  now  be  entered  to  obtain  actual  numbers. 
The  problem  formulation  is  usually  done  in  non-dimensional  coordinates 

X = jj,  t = ^ where  L and  a are  convenient  reference  length  and  sound 

speed,  respectively.  The  non-dimensional  time  interval  between  inter- 
section of  leading  and  trailing  waves  of  the  expansion  with  the  shock 
thus  yields  the  reference  length,  if  we  set  the  real  time  interval 
equal  to  the  pulse  duration.  The  numbers  as  yet  apply  to  a straight 
shock  tube,  chosen  to  give  the  desired  equivalent  downstream  conditions. 

In  converting  over  to  the  area-changed  shock  tube,  we  recognize 
that  the  driver  lengths  should  not  in  general  be  equal,  since  different 
pressure  ratios  are  involved  across  the  diaphragm.  Thus,  velocities 
and  wave  speeds  must  be  different;  and,  the  wave  speed  in  the  area- 
changed  tube  must  slow  down  in  the  steady,  subsonic  expansion  to  match 
the  downstream  conditions  set  for  the  equivalent  straight  shock  tube. 

The  wave  diagram  in  Figure  6 shows  a possible  wave-speed 
comparison  (as  seen  from  numbers  to  be  given  later)  and  illustrates 
how  a simple  conversion  of  tube  lengths  may  be  made:  we  demand 

that  the  real  times  be  equal,  for  both  straight  and  area-changed 
shock  tubes,  for  the  leading  wave  of  expansion  to  be  reflected  from 
breech  and  to  travel  back  to  the  diaphragm  station.  The  flows  to 
the  leading  waves  are  then  identical  beyond,  neglecting  the  time 
required  in  traversing  the  short  transition  piece  (which  implies  an 
assumption  on  the  length  of  this  piece). 
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This  particular  calculation  may  be  done  analytically  for  the 
intersection  point  of  leading  with  trailing  wave  of  the  reflected 
expansion8  and  wave  progress  then  continued  at  constant  speed  to 
the  diaphragm  station  X = 0 at  t . 

from  the  nondimens ional  times 


At 


a .At 
4 


the  real  times  for  the  arrival  at  X = 0 are  equated 


At 


) 

a4  / str 


At 


) 


c 


to  arrive  at  driver  length  conversion. 

C.  Sample  Calculated  Wave  Speeds  and  Results 


The  actual  test  conditions  are  used  in  the  calculations 
outlined.  As  a specific  example,  we  choose  the  test  condition 
p9 t = 73  kPa  (10.6  psi).  For  this  case  the  previously  obtained  flow 

parameters  are  given  to  an  equivalent  straight  shock  tube: 


M = 1.27 


Va  = 


7/5 


P2/Pl  = l«72 

u?  = 0.4022a^  = u^ 

M3  = u3/a3  =0.5  (M3a  = 0.6  for  A4/A  = 0.444) 


a3  = 0.8044a1 


With  these  conditions,  the  driver  sound  speed  is 

a.  = 0.88484a. 

4 1 


8.  R.  K.  Lobby  "On  the  Length  of  a Shook  Tube, " UTIA  Report  5, 
Institute  of  Aerophysios , University  of  Toronto , 1950. 
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Then  the  wave  speeds  for  the  reflected  expansion  are 


C 5 - u + a 
+ 5 5 

C+2  = u2  + a2 
C+y  = u?  + a? 
(C_4  h u4  - a4 
(C  3 = u,  - a. 


(0.4022  + 0.8044)a1  = 1.2066a 
0.72396  a, 

(0.4022  + 1.083)a1  = 1.4852a 
(-0.0610  + 0.9911^  = 0.9301 

- 0.88484  a ) 

- 0.40223^ 


With  these  results,  one  is  in  position  to  follow  the  rarefaction 
wave  from  diaphragm  break  to  its  interaction  with  the  primary  shock. 


A characteristics  calculation  with  desk  calculator  yields  the 
extent  of  the  breech  reflection,  for  leading  (£)  and  trailing  (t) 
waves : 


t = 1.1301 

al  1 

X = - 1.0,  (t  = f here) 

t = 2.15147  4 

The  intersection  point  of  reflected  leading  wave  with 
trailing  wave  is  found  to  be 

X3  = - 0.6094;  x3  = 1.5151. 

Continuing,  the  reflected  leading  and  trailing  waves  moving  with 
velocities  C+3  and  C+5,  respectively,  overtake  the  contact  surface 
K moving  with  velocity  u7.  We  have  bypassed  a detailed  (but  more 

accurate)  characteristics  calculation,  and  have  followed  both  waves' 
arrival  only  up  to  the  leading  wave's  intersection  with  K,  at  X : 

K 

tk  % = 3.0302 

XK  = 1.21876, 

tk  t = 5.21607 

This  procedure,  as  seen  by  inspection  of  the  wave  diagram,  Figure  5, 
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Should  give  ;i  lower  estimate  on  the  actual  trailing  wave's  inter- 
section point  with  K niul  hence  an  upper  estimate  to  he  made  on  the 
tube  lengths. 

Beyond  the  contact  surface  K,  the  transmitted  leading  and 
trailing  waves  of  the  rarefaction  move  to  intersection  with  the  shock, 
with  velocities  C 2 and  C 7 , respectively . The  procedure  for  this 

interaction  is  again  to  follow  the  rarefaction  waves  only  up  to  the 
leading  wave's  intersection  witli  shock.  Here  also  a detailed 

characteristics  calculation  is  bypassed  in  the  interest  of  simplicity. 
Similarly,  a lower  estimate  on  actual  intersection  position  and  time 
is  expected  for  the  trailing  wave.  The  results  then  are 

t . = 15.2495 
s , a 

X.  = 19.367, 

S t = 24 . 7 28 

s,  t 

The  (nondimensional)  interaction  time 


a^  At 


thus  fixes  the  desired  peaked  wave's  duration  (which  may  be 
transmitted  with  some  positive  pressure,  depending  on  conditions, 
nevertheless).  Choosing  the  wave  duration  At  = 50  ms  then  determines 
the  required  (straight  tube)  driver  length: 

a 50  (10-3) 

L = — 

4 str  At 


L str  = 1.77m  (5.80  ft.)(a  = 335  m/s  or 

1100  ft/s) 

The  test  wall  is  placed  at  X , where  the  leading  wave  of  the  rare- 
faction intersects  the  primary  shock  and  determines  the  driven  tube 
length  L^: 

Lj  = Xs  L4  = 19.367  (1.77m) 

L = 34.3m  (112.3  ft)  ; 

as  we  have  remarked  earlier,  as  determined  is  probably  an  upper 
estimate. 
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These  numbers  are  as  yet  based  on  the  straight  shock  tube  with 
downstream  conditions  identical  to  the  desired  conditions.  To 
maintain  this  identity  through  wave  speeds  and  interaction  points  in 
the  wave  diagram,  we  asked  that  the  real  times  be  equal  after 
diaphragm  break  - for  straight  and  for  area-changed  tube  - for  the 
leading  waves  of  the  expansions  to  travel  back  to  the  diaphragm 
station.  Thus,  using  the  numbers  appropriate  for  the  straight  tube 
and  for  the  area-changed  tube  for  the  case  considered  in  the 
beginning  of  this  section  and  following  the  procedure  outlined 
earlier,  we  arrive  at 


At 


L 


At 
c 

4c  a. 
4c 


At 


4 str  a. 


str 

str 


and  with  numbers  inserted 


or 


^ _ , 1.756  _ , 1.7747 

4c  aj  4 str  0.88484 


L.  = 1.142  L. 

4c  4 str 


For  L.  _ = 1.77m  (5.80  ft), 

4 str 

L = 2.02m  (6.62  ft) 

Then,  having  forced  this  partial  simultaneity  on  driver  wave  motions, 
we  are  assured  that  the  downstream  behavior  in  the  two  shock  tubes 
will  be  reasonably  matched.  Hence  the  driven  tube  length  is  as 
computed-  34.3m  (112  ft). 


IV.  CONCLUSIONS 

Using  an  extension  of  the  normal  straight  shock  tube  operating 
cycle,  we  have  determined  appropriate  tube  lengths  and  driver 
pressures  for  requisite  operation/modification  of  an  existing  BRL 
2.4m  (8')  shock  tube  to  permit  full  scale  testing.  Walls  and 
structures  of  3.7m  (12')  height  are  acceptable  for  simulated  blast 
loadings  of'vJO  kPa  (10  psi)  having  pulse  durations  of  50  ms.  Other 
loadings  and  the  appropriate  tube  lengths  are  derivable  without 
difficulty  using  the  procedures  developed. 

The  results  are  based  on  idealized  shock  tube  operation,  with 
approximations  introduced  to  obviate  the  need  for  an  extensive 
computer  program.  In  an  actual  design  to  the  full  sized  tube,  one 
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should  expect  departures  from  such  ideal  operation  due  to  boundary 
layer  growth  and  imperfect  diaphragm  effects  such  as  non-instantaneous 
opening  and  loss  of  energy  to  the  opening  process. 


Due  to  the  relatively  short  length  of  driver  required,  a much 
smaller  prototype  tube  with  required  area  change  ratio  could  be 
constructed  from  BRL  existing  shock  tubes,  and  the  wave  characteristics 
, determined  empirically  and  compared  with  calculations,  to  see  what 
modifications  should  be  incorporated  into  the  full  sized  facility. 
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Figure  1.  BRL  Dual  Shock  Tube  Facility 


Figure  2.  Wave  Diagram  of  Idealized  Shock  Tube  with  Area 
Change  at  Diaphragm 


26 


Figure  3.  Graphical  Solution  to  Equation  (8):  Curves  Z^a 
and  as  a Function  of  Mach  Number 


A4/A,=  0.444 


Figure  4.  Driver  Pressure  - Testing  Pressure  for  2.4m:  3.7m  ( 8 ' : 1 2 1 ) 

Area  Changed  Shock  Tube 
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